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Abstract. The set of quantum Gaussian channels acting on one bosonic mode can be 
^ ■ classified according to the action of the group of Gaussian unitaries. We look for bounds 

O" 1 on the classical capacity for channels belonging to such a classification. Lower bounds can 

be efficiently calculated by restricting to Gaussian encodings, for which we provide analytical 
CN| ■ expressions. 
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1. Introduction 

Gaussian processes are ubiquitous in physics, mathematics and information theory JH. In 
information theory, Gaussian channels are used to model noisy communication lines where 
the noise is described as a Gaussian process. A central problem in this field is to determine 
the capacity of the communication line, that is, the maximum rate at which information can 
be reliably transmitted via the noisy channel, with asymptotically vanishing probability of 
error [2]. The seminal work of Shannon, besides founding the theory of information, has 
provided the expression for the capacity of the Gaussian channel Q. On the other hand, 
the new research field of quantum information science has represented a fertile extension of 
information theory, with a heritage of open problems, challenges, and physical insights into 
the foundations and modern applications of quantum physics flUHl. 

Similarly to its classical counterpart, a quantum Gaussian channel is a mathematical 
model for a communication line associated with a noise process with Gaussian characteristic 
O (see also Q, where this type of maps were considered in view of dynamical invariants [H 
for quantum systems with quadratic Hamiltonians). Actually, differently from the classical 
case, there are several nonequivalent notions of capacities for a quantum channel, depending 
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on whether the aim of the communication protocol is to transfer classical or quantum 
information, and on whether entanglement is used to assist the transmission [5J. Physically, 
the model of quantum Gaussian channel can be applied in the context of continuous variable 
quantum systems, that is, to the case of bosonic systems, to model linear attenuation and 
amplification processes. Typical realizations include the quantum electromagnetic field and 
atomic ensembles. 

In the present contribution we consider the problem of evaluating the classical capacity, 
that is, the maximum rate of reliable transmission of classical information, through a quantum 
communication line modeled as a quantum Gaussian channel. This problem is formulated as 
an optimization problem, where a suitable entropic function — the Holevo information — has 
to be maximized over all possible protocols to encode classical information into an ensemble 
of quantum states [9J. The main technical difficulty is due to the fact that the Holevo function 
might be non additive ifTOll . which means that the optimization has to be carried over the whole 
unbounded set of encoding protocols, including those defined by ensembles of entangled 
quantum states. While the additivity of the Holevo function has been proved for a class of 
discrete quantum channels ifTTTl . nothing is known in general about continuous ones (to which 
the quantum Gaussian channels belong to). Then, the problem is dramatically simplified 
by restricting to the one-shot capacity, i.e., considering only ensembles of separable states. 
Actually it is further simplified if we also accept the conjecture that the optimal encoding 
ensemble for a Gaussian channel is solely made of Gaussian states [fl2|. Here, resorting on 
these assumptions, we compute lower bounds on the classical capacity of quantum Gaussian 
channels. 

The paper proceeds as follows. In Sec. [2] we briefly introduce the notion and the basic 
properties of quantum Gaussian channels; in Sec. [3] we considered the problem of estimating 
the capacity of Gaussian channels and compute lower bounds on the capacity of one-mode 
Gaussian channels; we end with conclusions in Sec.Hl 



2. Gaussian states and channels 



The most natural way to introduce the notion of quantum Gaussian channels is via that of 
Gaussian states. We consider the case of one-mode Gaussian channels, whose input and 
output quantum systems are described by a canonical pair R = (Q, P), obeying canonical 
commutation relations [Q,P] = i (here and in the following we set h — 1). A state of the 
quantum system is described by a density operator p, from which one writes the characteristic 
function 



Tr 



pW{Z) 



(1) 



where Z = (X, Y) is a two-dimensional real vector, and W(Z) = e is the Weyl operator. 
Gaussian states are, by definition, those with Gaussian characteristic function, i.e., 



exp ^Z Z T - izVZ T 



(2) 
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where 



Tr 



pR 



is the vector of first-moment of the canonical variables Q, P, and 

q2 qp+pq 



V= Tr 



P 



QP+PQ 
2 



2 
p2 



(3) 



(4) 



is the covariance matrix (CM). In addition to be positive definite, V > 0, the uncertainty 
relation imposes the following constraint on the CM associated with Gaussian states lfT3l : 



det(V)>-. 



(5) 



The latter condition can be generalized to the case of multimode Gaussian states, in terms of 
the symplectic invariants of the CM [13]. Actually, in the case of one-mode Gaussian states, 
the determinant of the CM is the only symplectic invariant [fl~4|l . In the following we will 
evaluate the capacity of Gaussian channels in terms of entropic functions. We hence recall 
that the von Neumann entropy, S[p) = — Tr(p log 2 p) (measured in qubits), of a Gaussian state 
is a function of the symplectic invariants of its CM. For the one-mode Gaussian states we have 

S[p] = h (VditOO) = S [V] , (6) 

where we have introduced the function h(x) = (x + 1/2) log 2 (x + 1/2) — (x — 

i/2)io g2 (x- 1/2) ma. 

Quantum Gaussian channels are those quantum dynamical maps [IUEE61 which transform 
Gaussian states into Gaussian states [6J. It follows that a Gaussian channel is identified by its 
action on the CM and on the vector of first moments. The preservation of the Gaussianity of 
the characteristic function implies that a Gaussian channel is identified by a triplet (d, T, N), 
whose action on the vector of first moment and the CM is as follows: 



Z ^ TZ + d , 
V h-> TVT T + N . 



(7) 
(8) 



where N is a 2 x 2 symmetric, positive semi-definite, matrix, T is a 2 x 2 real matrix, and d is 
a 2-component real vector. Moreover, the condition of complete positivity on the dynamical 
map is characterized by the inequality — stated in terms of the symplectic invariants of the 
matrices T, N — 



det(N) ^ 



det(T) - 1 



(9) 



Among the set of Gaussian channels, a remarkable role is played by Gaussian unitary 
transformations. These are characterized by triplets of the form (f , §, O), where f is a vector, 
§ is a symplectic matrix, and O denotes the null matrix. The composition, from the left and 
from the right, of the Gaussian channel with the Gaussian unitaries (f^, E>a, ©), (fs, §b, O), 
transforms the associated triplet according to 

d ^§ B (Tf A + d) + f B , (10) 

T^SbTS^, (11) 

N i — ^ § B N§X • (12) 
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B 2 


I 




{1} 


C (Att) 




(l-r)(n + 1/2)1 


(0,1) 


C (Amp) 




(r- l)(n+ 1/2)1 


(l,oo) 


V 




(1 -r)(n + 1/2)1 


(-oo,0) 



Table 1. Equivalence classes of one-mode Gaussian channels. Each class is identified by the 
matrices T c , N c . I and O respectively denote the identity and null matrices, Z = diag(l, —1). 

From the point of view of information theory, the composition from the left and from 
the right of a Gaussian channel with Gaussian unitaries corresponds to unitary pre-processing 
and post-processing of the channel. The Gaussian channels can be classified according to 
equivalence up to Gaussian unitary transformations: Two Gaussian channels are said to be 
equivalent if there exist Gaussian unitary pre-processing and post-processing mapping one 
channel into the other. For one-mode Gaussian channels, different equivalence classes can 
be identified [fTTl . which are summarized in Tabled! For each class, one can pick up a 
representative Gaussian channel associated with a triplet in the canonical form (0, T c , N c ), 
where = S B (Tf A + d) + f B , T c = § B T§ A , N c = § B N§J, and both matrices T c , N c are 
diagonal (see Table CD for details). Besides the rank-deficient classes A 2 , B%, each class is 
identified by the range of the parameter r := det(T). Within each class, a canonical form is 
characterized by the values of r and of the non-negative parameter n. The canonical form can 
be physically represented in terms of two-mode Gaussian unitary transformations, where an 
ancillary mode is introduced in a thermal state ifPTll : The class C describes linear attenuation 
(for r < 1) and amplification (for r > 1) processes, which are commonly modeled in terms 
of beam-splitters and linear amplifiers; The class V is the conjugate channel of the linear 
amplifier; Class Ai models the erasure process, and can be represented as the limit of the 
attenuating channel for r — >• 0. Class B 2 models the action of classical noise, its Stinespring 
dilation can be obtained with the help of two ancillary modes [?]. Finally, the classes A 2 , £>i 
present special features, they are characterized by rank deficient canonical matrices and can 
be obtained as singular limits of the previous cases. 

3. Classical capacity of one-mode Gaussian channels 

An encoding procedure is identified by a map x — >■ p x , which assigns a quantum state to each 
value of a stochastic variable, distributed according to a probability density p x . 

If the quantum states are subjected to the action of a quantum channel £, the receiver has 
to extract information from an ensemble composed of the states S(p x ). According to [9J, the 
maximum information, measured in bits, that can be extracted from the enseble is given by 
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the Holevo information: 



X = s 



dxp x £ (p x ) 



- / dxp^S 



(13) 



where S is the von Neumann entropy. 

The capacity of a quantum Gaussian channel £, associated with the triplet (d, T, N), is 
the supremum of the quantity (fT3l) over all possible encoding procedures. Among them, we 
consider Gaussian encodings, in which the value x of the stochastic variable is encoded in 
Gaussian states with CM V and first-moment d = x. Furthermore, we assume a Gaussian 
form for the distribution p x , with zero mean and CM M, which is symmetric and positive 
semi-definite. In these settings, we have that 



£(p x ) = S [TVT T + N] , 



S 



dxp x £(p x ) 



S [T(V + M)T T + N] 



(14) 



(15) 



However, the maximization of the Holevo information, even when restricted to Gaussian 
encodings, can lead to an infinite value of the channel capacity. This is due to the fact that the 
relevant Hilbert space is infinite-dimensional (and, hence, it can in principle store an infinite 
amount of classical information), and that the manifold of Gaussian states is not compact. An 
effective cutoff in the Hilbert space can be introduced by imposing a constraint on the average 
mean energy involved in the encoding process 0: 

^Tr(V + M)^£. (16) 

The optimization over Gaussian encodings furnishes a lower bound on the capacity of 
the Gaussian channel: 



C = max {S [T(V + M)T T +N] — S [TVT T + N] } 



(17) 



where the maximum is over the CMs V, M subjected to the uncertainty relation © and to the 
energy constrain (U~6l) . In the following the maximization problem will be solved by using the 
Karush-Kuhn-Tucker (KKT) method 0, which generalizes the Lagrange method to the case 
of constraints expressed by inequalities. 

Before proceeding with the evaluation of the lower bounds on the capacity for the 
Gaussian channels belonging to different equivalence classes, we remark that the von 
Neumann entropy is invariant under symplectic transformations, that is, 

S [fVT T + N] = S [§ B TVT T §J + § B N§£] , (18) 

for any symplectic matrix which is interpreted as a post-processing Gaussian unitary 
transformation. It follows that the function (fT7l) can be rewritten as follows 



C = max {S [T C S A (V + M)S^Tj + Nj - S [T c § A V§^Tj + N c ] } . (19) 



\ Clearly, other choices are possible, e.g., to constraint the output energy. However, since the Holevo 
information (13[ is obtained by optimizing over all the possible measurement on the channel output, we do 
not assume any constraint on the output of the channel. 
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This expression is written in terms of the canonical matrices and the pre-processing symplectic 
matrix E>a- The latter can be further decomposed according to Euler decomposition: 

§ A = RUM' , (20) 

where R, R' are orthogonal matrices and D = diag(r 1 / 2 , r -1 / 2 ). Since the orthogonal matrix 
R' preserves the constrain (fT6l) . it can be eliminated by redefining the CMs V, M, yielding 

C = max {S [T C RD(V + M)(T C RD) T + N c ] - S [T c ROV(T c R©) T + N c ] } .(21) 

3.1. Class A\: erasure channel 

Since T c = O, for this class of channels we trivially have 

C = max {S \N C ] - S [N c ] } = . (22) 

V,M 

Actually, by allowing non-Gaussian encoding we can similarly prove that the channel capacity 
itself vanishes. 

3.2. Class A2 

For the channels belonging to this class, T c = (I + Z)/2 = diag(l, 0) and N c = (n + 1/2)1. 
Since the matrix T C RD appearing in Eq. (|2TI) has rank one, its singular value decomposition 
reads 

T C RD = V^M'TcM" , (23) 

where W, W are orthogonal matrices, and y/t is the (non-vanishing) singular value of T C MD. 
Absorbing the orthogonal matrix M." into the definition of the CMs V and M, we obtain the 
following expression for the lower bound: 

C = max {S [tTJV + M)T C + N c ] - S [tT c VT c + N c ] } . (24) 

V,M L J 

By introducing the parameterization V = diag(s _1 /2, s/2), Eq. (|24l) can be rewritten as 
the maximum of a function of s: 



where 



C = max {h{u A2 ) - h{u M )} , (25) 

s£s_,s+ 



9 A2 = Mt (2E + 1 - -) +n+ -j \ n + - J , (26) 
V M = \ l*V + n + I J (n+l), (27) 



and s± = 2E + 1 ± ^(2/^ + l) 2 - 1. 
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3.3. Class Bi 

For the elements belonging to this class, T c = I and N c = (I + ^)/2 = diag(l, 0), Eq. (|2TI) 
reads 

C = max {S [MB(V + M)(R©) T + Nj - S [T c R©V(T c MO) T + N c ] } (28) 
= max { S I" V + M + ©" 1 M T N C M©" 1 1 - S [V + D _1 R t N c M t O _1 1 ) .(29) 

V,M L L J L J ' 

The matrix ID) _1 R T N C M]> _1 is symmetric with rank one, hence we can write it as follows: 

D -1 R T M C MD -1 = nR'N c R' T , (30) 

where n is its (non-vanishing) eigenvalue, and R' is orthogonal. Absorbing R' into the 
definition of the CMs V, M, we obtain 

C = max {S [V + M + nN c ] - S [V + nN c }} . (31) 

V,M 

Using the parameterization V = diag(s -1 /2, s/2), Eq. ([311 yields 

C = max {/i(P B J - h(u Bl )} , (32) 

where 

^ = ^(^y + n ) ( 2E + 1 -' 1 y) ' (33) 

^ = vz+^| ■ (34) 

5.4. C/aw C: attenuating and amplifying channels 

For this class we have T c = y/rl, N c = \ 1 — r| (n + 1/2)1, where r e (0, 1) for the attenuating 
channels, and r £ (1, oo) for the amplifying ones. Equation (|2TI) simplifies as follows 

C = max {S [rB(V + M)D + N c ] - S [rDVD + N c ]} . (35) 



First of all, we notice that the matrices D, N c in Eq. (1351) are all diagonal, hence the 
maximum under the constraint (TT6l) is reached in correspondence to diagonal CMs V, M. The 
constrained optimization yields the analytic solution 



C = h 



- h 



r(E + i) + |l-.|^ + ^ r + r " 

T 11 if 1 

- + 1 -r n + - 
2 1 1 V 2 



(36) 



The latter is obtained in the region of parameters for which 



± J 1 — n — (n+o)--s-^°- ( 37 ) 



2/ r 2 V 2/ 2 
Outside this region the KKT method does not lead to an analytical expression for 
the capacity lower bound, which can however still be evaluated numerically by solving 
a transcendental equation. The existence of qualitative differences in the solution of the 
optimization problems according to the range of the parameters has been put forward and 
analyzed in details in H9]| (see also |20l and ftZ\\ ). 
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3.5. Class £>2-' additive noise channels 

For the channels belonging to the equivalence class of the additive noise channel, T c = I, 
N c = nl, and Eq. $2A} reads 

C = max {S [B(V + M)D T + N c l - S [DVD T + N c l } , (38) 

— V,M L L J L J J 

This expression can be obtained from (1351 ) by replacing r with 1, and |1 — r\(n + 1/2) with 
n. Under the conditions 

E + ^±n— ->0, (39) 

the optimization yields the following analytical expression 

C = h^E + ^ + n r -±^j -h(n+^j . (40) 

Outside the region of parameters defined by Eq. (1391) the KKT method does not lead to a 
closed form for C_, see also f[T8l . 

3.6. Class T>: Conjugate of the amplifying channel 

For the channels belonging to this class, T c = ^—rZ, and N c = (1 — r)(n + 1/2)1, where 
r E (— oo, 0) and Z = diag(l, —1), yielding the following form for Eq. (|2T|) : 

C = max {S [|r|ZKB(V + M)(MD) T Z + N c ] - S [|r|ZMOV(MD) T Z + N c ] } .(41) 

Since the von Neumann entropy S is a function of the determinant of its argument, we can 
write 

C = max {S [|r|RB(V + M)(MD) T + ZN C Z] - S [|r|MDV(MD) T + ZN C Z] } 
= max {S [|r|D(V + M)D T + N c ] - S [|r|DVD T + N c ] } . (42) 



The latter expression is identical to Eq. (I35I) . upon replacing r with |r|, and |1 — t\ with 1 + \r\ 
Hence the solution of the maximization problem is obtained as for the class C, yielding 



C = h 
- h 

under the conditions 



r\[E + l)+(l + \r\)(n+ 1 -^ 



2J y 1 " V 2/ 2 

^ + (i + H)(n + i 



(43) 



^ 1\ (1 + Irhr-r- 1 / l\ r Tl 



2/ r 2 V 2/ 2 



4. Conclusions 



We have considered the whole family of one-mode quantum Gaussian channels, which can be 
classified into equivalence classes, and, by using Gaussian encoding, we have evaluated lower 
bounds on the classical capacity under input energy constraint. 
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For small enough values of the energy the optimal encodings can be evaluated 
numerically, while a closed form for the lower bound has been obtained when the allowed 
energy is above a certain threshold. The presence of two different regimes for the optimization 
problem has been studied in details in |fl9l (see also lfT8ll20ll2TT0 . The region of parameters 
allowing an analytical expression for the lower bound on the capacity includes the limit 
E — > oo. Using the fact that h(x) ~ log 2 a; + log 2 e for x — > oo, we have the following 
high-energy limits of the capacity lower bounds: 



The Holevo function (l2~TI) assumes different forms, depending on the equivalence class 
which the Gaussian channel belongs to. However, interestingly enough, the value of its 
maximum is a function of the channel itself and is not uniquely determined by the property of 
belonging to a given class. 

This is reminiscent of other problems in quantum information science for which the 
orbits of the action of the relevant unitary group do not provide a complete characterization 
of the theoretical framework (|22|. 
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